Dynamics of a Bright Soliton in Bose-Einstein condensates with Time-Dependent 
Atomic Scattering Length in an Expulsive Parabolic Potential 
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We present a family of exact solutions of one-dimensional nonlinear Schrodinger equation, which 
describe the dynamics of a bright soliton in Bose-Einstein condensates with the time-dependent 
interatomic interaction in an expulsive parabolic potential. Our results show that, under the safe 
range of parameters, the bright soliton can be compressed into very high local matter densities by 
increasing the absolute value of atomic scattering length, which can provide an experimental tool 
for investigating the range of validity of the one-dimensional Gross-Pitaevskii equation. We also 
find that the number of atoms in the bright soliton keeps dynamic stability: a time-periodic atomic 
exchange is formed between the bright soliton and the background. 
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With the experimental observation and theoretical 
studies of Bose-Einstein condensates (BECs)pJ, inten- 
sive interests have been paid to the nonlinear excitations 
of the atomic matter waves, such as dark[2j and bright 
solitons^, 0,0>@1- It is believed that the atomic matter 
bright solitons are of primary importance for developing 
concrete applications of BEC in the future, so it is of 
interest to develop a new technique which allows us con- 
struct a particular soliton with the assumed peak mat- 
ter density. One possibility is to vary the interatomic 
interaction by means of external magnetic fields. Re- 
cent experiments have demonstrated that the variation 
of the effective scattering length, even including its sign, 
can be achieved by utilizing the so-called Feshbach res- 
onance Q. This offers a good opportunity for manipu- 
lation of atomic matter waves and nonlinear excitations 
in BEC H, @ . In Ref. , it has been demonstrated 
that the variation of nonlinearity of the Gross-Pitaevskii 
(GP) equation via Feshbach resonance provides a power- 
ful tool for controlling the generation of bright and dark 
soliton trains starting from periodic waves. Besides, a si- 
nusoidal variation of the scattering length has also been 
used to form patterns such as Faraday waves [llj . or as 
a means to maintain BECs without an external trap in 
two dimensions |l2| . 

In this letter, we present a thorough analysis on the 
dynamics of a bright soliton of BEC with time-varying 
atomic scattering length in an expulsive parabolic poten- 
tial. Our study is greatly facilitated by the so-called Dar- 
boux transformation |l3|, by which we can directly con- 
struct the exact solutions of one-dimensional(lD) nonlin- 
ear Schrodinger equation(NLSE). Under the safe range of 
parameters, the bright soliton in BEC can be compressed 
into very high local matter densities by increasing the 
absolute value of atomic scattering length with Feshbach 
resonance. During the compression of the bright soliton 
in BEC, the number of atoms in the bright soliton keeps 



dynamic stability and the exchange of the atoms between 
the bright soliton and the background of density exists. 

At the mean-field level, the GP equation governs the 
evolution of the macroscopic wave function of a 3D BEC. 
In the physically important case of the cigar-shaped 
BEC, it is reasonable to reduce the GP equation into 
ID NLSEHElOil, 



M{x,t) d 2 ip{x,t) 



dt 



dx 2 



+ 2a(t)\^(x,t)\ 2 ^(x,t) 
+ -\ 2 x 2 i>{x,t) = 0. 
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In Eq. (JTJ, time t and coordinate x are measured 
in units 2/luj_ and ax, where a±_ = (h/muJx) 1 ^ 2 and 
o-o = (fr/muo) 1 / 2 are linear oscillator lengths in the trans- 
verse and cigar-axis directions, respectively. uj± and loq 
are corresponding harmonic oscillator frequencies, m is 
the atomic mass and A = 2|ti>o|/ w J_ <C 1. The Feshbach- 
managed nonlinear coefficient reads a (t) = \a s (t)\/aB = 
go exp (Xt)(aB is the Bohr radius) [la. Il7|. The normal- 
ized macroscopic wave function ip (x, t) is connected to 
the original order parameter ^>(r,t) as follows, 



¥(r,t) 



y/2iniBa 
x exp —iu±t 
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From the viewpoint of stability, 3D and ID equations 
are very different. For a true ID system, one does not 
expect the collapse of the system with increasing number 
of atomsp|- However, it happens that a realistic ID limit 
is not a true ID system, with the density of particles still 
increasing due to the strong restoring forces in the per- 
pendicular directions. To avoid the collapse of the bright 
soliton j we must restrict our study of BEC to the 
safe range of parameters, in which the system becomes 
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effectively ID, i. e. the energy of two body interactions 
is much less than the kinetic energy in the transverse 
direction^: e 2 = a ± /£ 2 ~ N\a s \/a < 1 (£ is the heal- 
ing length) . The bright soliton in BECs has been created 
with the parameters ofiVs; x 10 3 , ujj_ — 2tt x 700Hz and 
uiq = 2irx 7 'Hz, a,fi na i = —Aclb for 7 Li 4], which provides 
the safe range of parameters. With the same experimen- 
tal conditions in Ref. Q and a s (t = 0) = — 0.25as, we 
can calculate e 2 = a±/£ 2 ~ N\a s \/a = 9.5 x 10~ 3 < 1. 
Then, the scattering length is increased in the form of 
a (t) = goexp(Xt). After at least up to 50 dimension- 
less units of time, the absolute value of the atomic scat- 
tering length turns to |a s (t)| = O.Sas, corresponding to 
e 2 = a±/£_ 2 ~ N\a s \/a = 3xl(T 2 < 1. Under the above 
conditions, the system is effectively ID. So the safe range 
of parameters can be described as follows: (1) with the 
same experimental conditions in Ref. 0; (2) ramp up 
the scattering length in the form of a (t) = go exp (Xt) 
within 50 dimcnsionless units of time. We also have to 
specify the terminology long-time dynamics. A unity of 
time, At = 1, in the dimensionless variables corresponds 
to 2/ui± real seconds. This means, for example, that 
for a BEC in a trap with transversal size of the order 
of a± ~ 1.5/im, a unity of the dimensionless time corre- 
sponds to 5.0 x 10~ 4 s. The lifetime of a BEC in today's 
experiments is of the order of Is, which is about 200 in 
our dimensionless units. 

The so-called 'seed' solution of Eq. Q can be chosen 
as follows, 



ipo (%, t) = A c exp 



where p c = koxexj>(Xt) 



Xt 



(3) 



(2 g0 ^-fcg)(cxp(2At)-l) 



4 ~ 2A 

and A c and fco are the arbitrary real constants. Wc 
perform the Darboux transformation [lflj ipi = V'o + 

( C+ f^ 102 exp (-At/2 - iXx 2 /4) to obtain the new 
solution of Eq. Q by taking Eq. © as the seed. Then 
we obtain the exact solution of Eq. (JJJ as follows: 



where 



A r +A. 



x exp 



(7 cosh 9 + cos (f) +i (a sinh 9 + (3 sin <p) 



cosh 9 + 7 cos ip 
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[(fco + fc.) A R - V^^sA/] [exp (2Xt) - 1] 



2A 



f = 



+Aflxexp (Xt) , 

[(fc + k s ) Aj + ^A S A R ] [exp (2Xt) - 1] 
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+ A/a; exp (Xt) , 



(5) 



and 



g A c (fc - k s + A/) 



(3=1- 
1 



2g Q A 2 c 
A ' 

^fgo~A c (A R - y/g^As 
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with 



A 



A = g A 



(A R 



1 (k s - fc )] 2 - AgoA 2 c = A R + iAi, 



9oA 2 s 



(k s - fc + A/) 
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where fc s is the arbitrary real constant. On the one hand, 
when A c = ko = 0, Eq. Q reduces to the well-known 
soliton solution: ip s = A s sech6* s exp(At/2 + iip s ), where 
$s = — y/go exp(Xt)A s x + y/go~k s A s [exp (2At) — 1] /A and 
Ps = fc — go A 2 [exp (2Xt) — 1] /2A. On the other hand, 
when the amplitude of this soliton vanishes (A s = 0) , 
Eq. J3J reduces to Eq. J2Jl. Thus, Eq. Q represents a 
bright soliton embedded in the background. Considering 
the dynamics of the bright soliton on the background, the 
length 2L of the spatial background must be very large 
compared to the scale of the soliton. In the real experi- 
ment , the length of background of BEC can be reached 
at least 2L =370/im. At the same time, in Fig. 1, the 
width of the bright soliton is about 21 = 2 x 1.4 /im=2.8 
\xm( a unity of coordinate, Ax = 1 in the dimensionless 
variables, corresponds to a±_ = (Ti/moj^) 1 / 2 = 1.4/im). 
So, we indeed have I -C L, a necessary condition for re- 
alizing our soliton in experiment. 

By utilizing the properties of Eq. JI) , we demonstrate 
that the manipulation of the scattering length can be 
used to compre ss a bright soiton of BEC into an assumed 
peak matter density. It has been reported that an abrupt 
change of the scattering length can lead to the splitting 
of the soliton with generating the new solitons. The frag- 
mentation of the soltion obviously decreases the numbers 
of atoms of the original soliton, which is undesirable for 
application^. However, in Eq. (gj, the change of the 
scattering length preserves the soliton from splitting new 
ones. For simplicity we assume fco = k s in Eq. and 
only consider the case of A 2 > AA 2 , for in the case of 
A 2 < AA 2 , a small perturbation for Eq. (@J may lead to 
the modulation instability [2l| • On the conditions above, 
Eq. (0} can be deduced into the following form: 



-A c + 6 2 



82 cos p — iA s sin p 



A s cosh 9 — 2 A c cos p> 



Xt 

x exp \ — + ip c , ■ 



(8) 



where 



g 8 2 xexp(Xt) 



gok S 2 [exp(2At) - 1] 



A 



g A s S 2 [exp(2Ai) - 1] 



2A 



y/Aj-~4Aj. 



(9) 
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For a better understanding, we plot Eq. 



in Fig. 
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Figure 1: The dynamics of Feshbach resonance managed soli- 
ton in the expulsive parabolic potential given by Eq. © 
The parameters are given as follows:A = 2 x 10 
A c = 1,A S = 2.4, fc = 0.03, 
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1, which shows the dynamics of the Feshbach resonance 
managed bright soliton in the expulsive parabolic poten- 
tial. As one can see from Fig. 1, with the increasing 
absolute value of the scattering length, the bright soliton 
has an increase in the peaking value, and a compression 
in the width. As a result, we can get a bright soliton 
with the assumed peak matter density. It is interesting 
to observe that in the expulsive parabolic potential, the 
bright soliton is set into motion and propagates in the 
longitudinal direction, instead of oscillating in an attrac- 
tive parabolic potential. The possibility of compressing 
BEC solitons into an assumed peak matter density could 
provide an experimental tool for investigating the range 
of validity of the ID GP equation. Since the quasi-lD GP 
equation only applies for low densities, it would indeed 
be interesting to see how far one can compress a soliton 
in a real experiment by increasing the a bsolute value of 
the scattering length. 

Inspired by two experiments 0, Q, we can design the 
compression of a bright soliton in BECs with the fol- 
lowing steps: (i) create a bright soliton in BECs with 



xlfp 



2tt x 700Hz and 



the parameters of N 
u>o = 2ir x 7 Hz, a s = — 0.25as for 7 Li. The main ef- 
fect of this expulsive term is that the center of the BEC 
accelerates along the longitudinal direction, (ii)under 
the safe range of parameters discussed above, ramp up 
the absolute value of the scattering length according 
to a(t) — goexp(Xt) due to Feshbach resonance, where 
A = 2|o;o|/w_l = 2 x 10~ 2 is a very small value. A unity of 
time, At = 1 , in the dimensionless variables corresponds 
to 2/ui± = 4.5 x 10~ 4 real seconds, (iii) after at least up to 
50 dimensionless units of time, the absolute value of the 
atomic scattering length turns to |a s (i)| = 0.8as, which 



is less than |o/j na j| = 4as. This means that during the 
process of compressing the bright soliton, the stability 
of soliton and the validity of ID approximation can be 
kept as displayed in Fig. 1. Therefore, the phenomena 
discussed in this letter should be observable within the 
current experimental capability. 

Furthermore, based on Eq. (JSJ, we find that when 
sinh.0 = , the peak matter density of bright soliton 
arrives at the maximum 



\i>\ 2 =e*p(\t) [A 2 C + 



5,A S 



A s — 2A C cos tp 



(10) 



and when cosh( 



the peak matter 



A c cos (p A s 

density of bright soliton arrives at the minimum 



A 2 C 5\ cos 2 ip 



|V>| =exp(At) A 



A 2 - AA 2 



cos 2 ip 



(11) 



This means that the bright soliton can only be squeezed 
into the assumed peak matter density between the min- 
imum and maximum values. In order to investigate the 
stability of the bright soliton against the variation of the 
scattering length in the expulsive parabolic potential, we 
obtain 



lim 



+ L 



-L 



\ij{x,t)\ 2 -\^{±L,t)\' 



dx 



28 2 



(12) 



which is the exact number of the atoms in the bright soli- 
ton against the background described by Eq. (JHJ). This 
indicates that during the process of the compression of 
the bright soliton, the number of atoms in the bright 
soliton keeps invariant. In contrast, the quantity 



lim 

L — >oo 
26 2 



\ip(x,t) -ip(±L,t)\ dx 



-L 



/.9o 



(1 + A C M cos ip) . 



(13) 



where 



4 arctan 



M = 



■J A„+2A C cos ip 
cos ip 



^JA 2 S - 4A 2 cos 2 ip 



(14) 



counts number of atoms in both the bright soliton and 
background under the condition of ip(±L,t) 7^ 0. Eq. 
(|13fl displays that a time-periodic atomic exchange is 
formed between the bright soliton and the background. 
As shown in Fig. 2(a), in the case of zero-background, 
i. e. A c = 0, there will not be the exchange of atoms, 
However, in the case of nonzero-background as shown 
in the Figs. 2 (b) and (c), the exchange of atoms be- 
tween the bright soliton and the background becomes 
more quickly with increasing the absolute value of the 
scattering length. The conclusion can be made that the 
number of atoms in the bright soliton in BEC keeps dy- 
namic stability against the variation of the scattering 
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Figure 2: Time-periodic atomic exchange between the bright 
soliton and the background given by Eq. (T^). The range of 
time is: (a) and (b) t=[0, 5], (c) t=[20, 25]. The parameters 
are given as follows: A = 0.02, go = 1, A s = 4.8, (a) A c = , 
(b) and (c) A c = 2.3. 



length in the expulsive parabolic potential. The consid- 
eration provided above implies that we should construct 
the bright soliton in BEC on the background of density 
described by Eq. ©. A question arises about the pos- 
sibility of creation of such a state experimentally. We 
notice that Eq. © will take the particular form at the 



time to = 2j m 



lows 



A S S 2 go 



1 



At 



, n = — 1, —2, —3. ..as fol- 



V> (a;, t)= — Ac exp I ^ J exp (iip c ) 

—182 exp sechf? exp (iip c ) , 

= - ^0 (X, t ) - Ipsoliton {%, to) ■ (15) 

which is the linear combination of Eq. Q and a bright 
soliton. So Eq. I|15[l means that Eq. (JSJ) can be generated 
by coherently adding a bright soliton into the background 
of density described by Eq. @ . 

In conclusion, we present a family of exact solutions of 
the nonlinear Schrodinger equation with the time- varying 
nonlinear coefficient in the expulsive parabolic potential. 
Our results describe the dynamics of Feshbach resonance 
managed bright soliton of BEC in an expulsive parabolic 
potential. Furthermore, under the safe range of param- 
eters, it is possible to squeeze a bright soliton of BEC 
into the assumed peak matter density, which can pro- 
vide an experimental tool for investigating the range of 
validity of the ID GP equation. We also find that the 
number of atoms in the bright soliton keeps dynamic sta- 
bility: the exchange of atoms between the bright soliton 
and the background becomes more quickly with increas- 
ing the absolute value of the scattering length. Recent 



developments of controlling the scattering length in the 
experiments allow for the experimental investigation of 
our prediction in the future. 
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